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Abstract 

Measurements by the Belle Collaboration of the cross section for inclusive J /ip production in 
e^e^ annihilation have been a serious challenge to current heavy-quarkonium theory. Especially, 
the measured cross sections for exclusive J/ip -\-r]c and inclusive J / ip -\- cc + X differ from nonrela- 
tivistic QCD predictions by an order of magnitude. In order to check if other available alternative 
theory can resolve such a large discrepancy, we calculate the cross section for inclusive J/tp + cc+X 
based on the color-evaporation model. As a phenomenological model, the color-evaporation model 
is still employed to predict cross sections for inclusive quarkonium production in various processes. 
Our results show that the color-evaporation-model prediction is even smaller than the nonrelativis- 
tic QCD prediction by an order of magnitude. The resultant color-evaporation-model prediction 
for J/^p + cc + X fraction in the inclusive J/ip production cross section is 0.049, while the empirical 
value measured by the Belle Collaboration is 0.82. 

PACS numbers: 13.66.Bc, 12.38.Bx, 14.40. Gx 
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Spin-triplet S'-wave charmonium states such as J/ip have been clean probes of both 
perturbative and nonperturbative feature of quantum chromo dynamics (QCD). Systematic 
studies based on first principles became possible after the introduction of the nonrelativistic 
QCD(NRQCD) factorization formalism 0|, an effective field theory of QCD. However, still 
there are several challenges to quarkonium physics research. One of the most interesting 
problems is regarding inclusive J/ip production in e^e~ annihilation at B factories. The 
Belle Collaboration has measured the cross section for exclusive J/ip + Vc production by 
observing a peak in the momentum spectrum of inclusive J/ip signals that corresponds to 
the 2-body final state J/ip + ric 2]. The measured cross section is by an order of magnitude 
larger than the predictions of NRQCD 0, 0|. Interesting proposals 0, 01 to resolve 
the problem were disfavored by a recent analysis by the Belle Collaboration jq. Having 
failed in explaining the result using lowest-order perturbative calculations, one may guess 
that higher-order corrections in strong coupling may be very large If it is true, 
perturbative expansion is not a proper method to predict the cross section. Recently, we 
proposed that the measurement of the cross section for the e~^e~ annihilation into four 
charm hadrons may provide us with a strong constraint in determining the origin of the 
large discrepancy jlO| and corresponding experimental analysis is being carried out. In the 
Belle analysis |2| they also reported the cross section for the inclusive J/ip + cc production at 
the center-of-momentum(c.m.) energy ^/s = 10.6 GeV. From the measured inclusive cross 
sections a{e+e- J/i) + D*+ + X) = 0.53 to^^ ± 0.14 pb and cr(e+e- ^ JM + + X) = 
0.87 

-a28 =1= 0-20 pb they extracted the cross section for J/ip + cc + XW^ based on the 
Lund model for fragmentation of a charm quark into the D mesons The resulting 

cross section is a{e'^e~ — > JZ- ^ + cc + X) = 0.87^^11 ± 0.17 pb 2^, which is larger than 
NRQCD predictions [ll by an order of magnitude. For example, an NRQCD 

prediction given in Ref. [1^ is about 0.07 pb. The Belle Collaboration also calculated the 
ratio a(e+e- J/V' + cc + X)/a{e+e- J/ip + X) = 0.82 ± 0.15 ± 0.14 0, This 
remarkable result reveals the fact that e^e~ —>■ J/ip + cc is the dominant source for inclusive 
J/ip production in e"'"e~ annihilation at B factories. Again, the result is remarkably 
larger than available NRQCD predictions. Two-photon mediated process e"'"e~ —>■ J / ip+cc of 
order has been calculated using NRQCD in order to check if the photon-fragmentation 

effect is significant like the process e+e- J/^ + J/iP 0,i. Unlike the exclusive double 
J/ ip production the photon- fragmentation contribution was found to be negligible compared 
to the single-photon process of order a^a^ in e+e- ^ J/V' + cc HI. Different approaches 
were also introduced using a large i^-factor and the nonperturbative quark-gluon-string 
model 19l |. 



In this paper, we calculate the cross section for inclusive e+e J/iIj + cc + X process at 



the cm. energy ^/s = 10.6 GeV using the color-evaporation model(CEM) |20l. l2lL l22l l23l|. 
Predictions of the quark-hadron-duality model, which is similar to the CEM, are avail- 
able \2^, 2^ and they underestimate the cross section for e+e" J/ip + cc like that of 
NRQCD. The quark-hadron-duahty model assumes that only a color-singlet cc state evolves 
into a J/ip, while the CEM sums over all possible spin and color states. Thus the CEM 
calculation will probe the importance of the missing color-octet component. The hard pro- 
cess for the inclusive production of J/ip + cc is approximated by e+e" — *■ cccc in leading 
order in strong coupling a^. Invariant mass of the cc pair evolving into the final J/ip is 
restricted from 2mc to open charm threshold 2m£). Multiplying the CEM parameter Fj/^ to 
the hard-scattering cross section, the CEM prediction for the inclusive J/ip + cc production 
cross section is obtained. The CEM parameter Fj/^ represents the probability of the cc pair 
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FIG. 1: Two topologically distinct Feynman diagrams for e {ki)e^{k2) c(pi)c(p2)c(p3)c(p4) 
among 8 diagrams. The pair c{pi)c{p2) is to evolve into the final J/V'- 

evolving into a J/ip. Our results show that the GEM prediction is significantly smaller than 
that of NRQCD, which already underestimates the empirical cross section by an order of 
magnitude. The process is well distinguished from many other hadronic processes, where 
GEM predictions are not that far away from NRQCD ones. Comprehensive reviews on the 
CEM can be found in Refs. pl,!!^. 

In the CEM the cross section for inclusive J/ip production is obtained by integrating 
differential cross section for cc + X over the invariant mass rricc of the cc pair from 2mc to 
DD threshold, 2mD- 

a{J/i^ + X)=Fj,.J dm,,, (1) 

where a,, is the cross section for e~^e~ cc + X. The CEM parameter Fj/^ in Eq. 
represents the probability for the cc pair to evolve into the final J/ip inclusively. This 
phenomenological parameter is fitted to the data from hadronic experiments which have a 
large number of data samples. Thus the uncertainties of the parameter come from parton 
distribution function and renormalization/factorization scale as well as the charm-quark 
mass rric. 

Before we proceed to consider the inclusive J/ip + cc + X production, we must address the 
point that the CEM is a model valid only for inclusive production of a heavy quarkonium. 
The parton cross section a,, in Eq. (Q) includes all possible spin and color states for a cc pair. 
The final quarkonium state is always accompanied by light hadrons because the quantum 
numbers for cc and the quarkonium can be different in general. Therefore, it is impossible to 
apply the CEM to predict the exclusive two-charmonium production in e"'"e~ annihilation. 

In order to apply the CEM to calculate the cross section for the inclusive J/ip + cc 
production in e"''e~ annihilation, we have to calculate the invariant-mass distribution of a 
cc pair created in the process e"'"e~ —>■ cccc + X and follow the way shown in Eq. ([Q). In 
leading order in strong coupling a^, cccc can be produced at order a'^a^. There are two 
topologically distinct Feynman diagrams generating two pairs of cc, which are shown as A^i 
and A^2 in Fig. ^a) and^b), respectively. Momenta for the involving particles are assigned 
as e~{ki)e~^{k2) — * c{pi)c{p2)c{p3)c{p4^) . The amplitude for the two diagrams shown in Fig.[T] 
are 

where s = {ki + A;2)^, Cc = | is the fractional electric charge of the charm quark, and a is 
the SU(3) color index for the virtual gluon. The vector indices a and (3 are for the virtual 
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photon and the gluon, respectively. We suppress the spin and color indices of the charm 
quarks in Eq. For i = 1 or 2 the tensors Hf^ in Eq. which are matrices in spinor 
space, are defined by 

Hf = /A(pi+p2+P3)7", (3a) 
Hf = 7"A(-p2-p3-P4)7^ (3b) 

where A{p) = {^+mc)/{p'^ — ml). There are 6 more Feynman diagrams that can be obtained 
from the two amplitudes Aii and by exchanging two charm quarks and two antiquarks, 
respectively, as 

Ms = -Pi^sMi, = -Pi^3-M2, 

-Ms = -P2^4-Mi, Me = -P2^4-M2, (4) 

Mr = +Pl^3P2^4Mu Ms = +Pl^3P2^4M2. 

where Pi^j is the operator exchanging two particles with momentum indices pi and pj shown 
in Fig. ^ The signs of Ms through A^g in Eq. are determined by the antisymmetry of 
Fermi statistics in exchanging identical fermions among the final-state particles. 

In the process e~^e~ cccc, there are four ways to pick a pair of cc, which will evolve into 
the final J/ip inclusively. We choose the cc pair with momenta pi and p2- Then we integrate 
over the invariant mass = a/ {pi + ^2)^ of the pair from 2mc to 2m£). The invariant 
mass = a/ (ps + Pa^ of the other pair is integrated over the whole phase space, which 
runs from 2mc to v^ — 2mi2. Because we distinguish two charm quarks and two anti-charm 
quarks simultaneously, we do not multiply the statistical factor (1/2)^ for identical particles 
in the phase space integration. Corresponding CEM prediction for the inclusive process 
e'^e^ cccc + X is calculated using the formula 

/■^'"■^ f da \ 

a{J/i) + cc + X) = Fj/^ / rfmi2, (5) 

J2mc \"'^12/ 

where the differential cross section in Eq. ^ is defined by 

r/rr 1 rVs-mu r 

i^rw^^L *"«/<'"''":^''«3Jp|lp;llP^El>'l^ (f) 

where M = X]i=i-Mi- The summation notation ^ in Eq. (jH)) stands for averaging over 
initial spin states and summation over final color and spin states. The solid angle dQ is 
for the three momentum P of pi + p2 in the e"'"e~ cm. frame. Remaining two solid-angle 
elements of the form dQ*j are for the three momentum p* in the rest frame of pi +Pj, where 
{i,j} = {1,2} or {3,4}. Integrating over the eight variables in Eqs. (0) and we get 
the CEM prediction for the cross section cr(e+e^ J/ip + cc + X). The cross section for 
e"'"e~ — > cccc can be obtained if we remove the Fj/^ in Eq. (0), replace the cut 2mD in Eq. © 
by — 2mc, and multiply the statistical factor (1/2)^ to consider two pairs of identical 
particles in the final state. The process e^e~ — > cccc has been considered in Ref. We 
compute the in Eq. (jH)) using REDUCE [2^ and carry out the phase-space integral 

in Eqs. (0) and (0) making use of the adaptive Monte Carlo routine VEGAS j29[. As a check, 
we carry out the same calculation using CompHEP j^^. Our analytic result for and 
numerical values for the total cross section agree with those obtained by using CompHEP. 

Let us summarize our theoretical inputs. Following recent quarkonium calculations in 
Refs. IE 0) we use the electromagnetic coupling a = 1/137 and the next-to- leading 
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order pole mass rric = 1.4 ± 0.2 GeV for the charm-quark mass. The strong couphng 
constant depends on the typical scale in the hard-scattering cross section. If we use the 
renormalization scale fi = 2mc with rric = 1.4 GeV, the value is as{2mc) ~ 0.26. The value 
is appropriate for the gluon-fragmenation contribution. However, a cc pair from different 
quark lines can evolve into the J/ip. In this case, the momentum transfer carried by the 
internal gluon can be as large as ^/i/2 = 5.3 GeV to get asU/s/2) ^ 0.21. The latter value 
has been used for the exclusive J/ip + r]c production in Ref . (3| because the momentum scale 
is exactly \/s/2. In our numerical analysis we take into account the uncertainty in the scale 
and allow the variation = 0.23 ±0.03. Because the cross section is proportional to a^, the 
uncertainty from the strong coupling is about ±30%. Finally, we have to fix the numerical 
value for the GEM parameter Fj/^ for directly produced J/ip. The GEM parameter i^}/^ 
for inclusive J/ ip production, which includes feed-down from higher charmonium resonances, 
has been fitted to the data from pp and pA collisions. They have small statistical errors 
but depend on patron distribution functions, renormalization /factorization scale, and m, 
A recent set of the values can be found in Refs. S H H Q . In Table 3.7 of Ref. ^ 



several values are given depending on parton distribution function and m^. The values 
are F}^^ = 0.0248 for = 1.4 GeV, 0.0229 for = 1.3 GeV, and 0.0144 ~ 0.0155 



for rric = 1-2 GeV. The fraction of direct J/^p is 0.62 as shown in Table 3.6 of Ref. |2'i 
The values for directly produced J/ip are then Fj/^ = 0.015 for rric = 1-4 GeV, 0.014 
for rric = 1-3 GeV, and about 0.01 for rric = 1-2 GeV. In Ref. j^, authors extracted the 
parameter using the open charm production in the GEM framework to estimate the QGD 
correction to the hadropro duct ion of charmonium. The values in Ref. ^] are consistent 
with the photoproduction data jl^, ll^. For rric = 1-45 GeV, the authors of Ref. jl^ Is^l 
found FJj"^ = (0.43 ~ 0.5)/9 which results in Fj/^ = 0.030 ~ 0.034. Roughly taking into 
account the difference between the two fits, we take Fj/^ = 0.025 ±0.010 for rric = 1-4 GeV. 
The uncertainty from the GEM parameter is about ±40%. The fits in Refs. Isp . l3^ favor 
larger values, while the fits in Ref. prefer smaller ones. As shown in Ref. |27l| numerical 
value for Fj/^ depends on rric. We take into account the rric dependence of the Fj/^ described 
above by making a linear fit to the values given in the Table 3.7 of Ref. [l^]. Resultant fit 
is Fj/^ = 0.025 ± 0.047 - 1.4). 

In Fig.|21we show the differential cross section da/dmi2 in Eq. (jU)) multiplied by the GEM 
parameter Fj/^ = 0.025 with = 1.4 GeV and = 0.23. The area of the shadowed region 
in Fig. 121 corresponds to the inclusive cross section cr(e"*"e~ J/ip + cc + X). If we increase 
the rUc, the cross section decreases because the left end point of the phase space for mi2 shifts 
toward the cut 2mD = 2 x 1.87 GeV and the height of the peak decreases. The cross section 
gets larger if we choose a smaller rUc. The rUc dependence of our GEM predictions for the 
inclusive cross section cr(e+e~ ^ J / ip + cc + X) at y/s = 10.6 GeV is shown in Fig. El The 
band represents the uncertainty ~ ±30% coming from the strong coupling as = 0.23 ±0.03. 
In Fig. El we use the mc-dependent GEM parameter Fj/^ = 0.025 ± 0.047 - 1.4) . 
Gonsidering roughly ±40% uncertainty in Fj/^ and ±30% uncertainty in strong-coupling 
dependence, our prediction has roughly ±50% uncertainty neglecting rric dependence, which 
is well described in Fig. 01 

In previous NRQGD calculations for e~^e~ J/ip + cc + X authors used larger values 
for the strong coupling as [13, EH For example, as{2mc) = 0.28 was used in Ref. (lj| 
to get the NRQGD prediction for the cross section cr(e+e" J /ip + cc + X) = 0.07 pb 
at ^/s = 10.6 GeV. In order to have a fair comparison with the NRQGD predictions let 
us replace our with = 0.28. Then our predictions given in Figs. El and El should be 
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FIG. 2: Differential cross section Fj/^da/dnicc in fb/GeV with respect to the invariant mass 
nice = of cc for e"*"e~ annihilation into cccc, where rric = 1.4 GeV, a = 1/137, = 0.23, 



and F 



0.025. The area of the shadowed region with 2mc < nice ^ 2m/) stands for the GEM 



prediction for (j(e+e — > J/V' + cc + X). 
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FIG. 3: The GEM prediction for the inclusive cross section a{e^e J/ip + cc + X) in fb at ^/s = 
10.6 GeV as a function of mc. The mc-dependent GEM parameter Fj/^ = 0.025+0.047 - 1.4) 
is used and the band represents the uncertainty from the strong coupling = 0.23 ib 0.03 only. 



increased by a factor of (0.28/0.23)^ ~ 1.5. The upper bound of our CEM cross section is 
about 12 fb which is less than the NRQCD prediction by a factor of 6. The underestimation 
of the cross section in the CEM may be due to the limited phase space for the invariant 
mass nice because of the small center-of-momentum energy. Further studies of invariant- mass 
spectra for hadronic processes will reveal if our guess is valid. 

Another important feature of our result is that the CEM prediction is less than that 
of the quark-hadron-duality model. The prediction of the quark-hadron-duality model is 
about 60 fb HI 113 which is similar to the NRQCD predictions HI El Q. The quark- 
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hadron-duality model collects only color-singlet contributions while the CEM includes all 
the color degrees of freedom. In the CEM the invariant mass of the cc pair is integrated 
over 2mc < rricc < '^frto-, which is less than the region 2mc < rricc < ^tud + (0.5 ~ 1 GeV) 
for the quark-hadron-duality model. If we increase the upper limit for the m^c by 1 GeV, 
the CEM cross section at = 1.4 GeV becomes 12 fb, which is still smaller than the 
quark-hadron-duality-model prediction by a factor of 5. A reason for this difference may be 
due to the long-distance factor Fj/^. We use the value for the Fj/^ extracted from hadronic 
collisions, while the quark-hadron-duality model uses a different method. 

Before closing our discussion, let us comment on the ratio a{e~^e^ J/ip + cc + 
X)/a{e+e- Jj^ + X) = 0.82 ± 0.15 ± 0.14 measured by the Belle Collaboration [l^. 
The empirical value is remarkably larger than the prediction of NRQCD. Unlike NRQCD, 
the nonperturbative factor Fjj^ in the CEM prediction for the ratio cancels and the ratio 
is depending only on and vrtf.. Therefore CEM prediction for the ratio should be more 
reliable than the absolute value for the cross sections. The inclusive Jjip production rate 
can be calculated by considering the parton process e"^e~ — > ccg in the leading order in a^. 
Note that e"'"e~ — >• cc channel is forbidden because rricc = 3> 2mD. According to our 
recent analysis js^, CEM prediction for the cross section is 

(T(e+e- ^ Jji) + X) ^ 95 fb, (7) 

where we use Fj/^ = 0.025 and as(2mc) = 0.26 with rric = 1.4 GeV. Resulting ratio becomes 

cr(e+e" J/%1) + CC + X) 



^ 0.062. (8) 

CEM 



a{e+e- J/ip + X) 

The empirical rate measured by the Belle Collaboration was obtained with the cut |p}/^| > 
2.0 GeV of the three- momentum of the J/ip in the e"'"e~ cm. frame. If we impose the 
same cut to our prediction, our prediction for the ratio becomes ~ 0.049, which is less than 
the value given in Eq. (jH)) by about 21%. The ratio is again significantly smaller than the 
experimental value ~ 0.82. 

In summary, we have calculated the CEM prediction for the inclusive cross section for 
e"'"e~ J/ip + cc + X u sing the CEM parameter Fj/^ fitted to the data from pp and 
pA collisions ji^ li^, |3l[ and from photoproduction [s^. Our CEM prediction for 
the cross section a{e^e~ J /ip + cc + X) is about 9.7 fb at = 1.2 GeV, 5.8 fb at 
rric = 1.4 GeV, and 2.2 fb at rric = 1.6 GeV. The cross section is, at least by a factor of 6, 
smaller than the NRQCD prediction, which is already smaller than the empirical value by an 
order of magnitude. Thus the color-evaporation model severely underestimates the inclusive 
J/ip + cc + X cross section in e"'"e~ annihilation measured by the Belle Collaboration. 
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